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The efficacy of isotope thermometry: examining in the S'-matrix approach 
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Isotope thermometry, widely used to measure the temperature of a hot nuclear system formed in 
energetic nuclear collisions, is examined in the light of S'-matrix approach to the nuclear equation 
of state of disassembled nuclear matter. Scattering between produced light fragment pairs, hitherto 
neglected, is seen to have an important bearing on the extraction of system temperature and volume 
at freeze-out from isotope thermometry. Taking due care of the scattering effects and decay of 
the primary fragments, a more reliable way to extract the nuclear thermodynamic parameters is 
suggested by exploiting least-squares fit to the observed fragment multiplicities. 
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Theoretical investigations on the equation of state 
(EOS) of infinite 0, 0, i, and finite @, @, 0, i, ©, E3] 
nuclear matter predict the existence of a liquid-gas (LG) 
type phase transition in these systems. This transition 
is thought to play an important role in nucleosynthesis 
in supernova explosion [111, 0] ■ Laboratory experiments 
in collisions between energetic nuclei appear to reveal 
signals of LG phase transition in hot finite nuclear sys- 
tems m q m. 

Proper identification of such a tran- 
sition, however, depends on the reliable measurement of 
the thermodynamic observablcs. In particular, temper- 
ature plays a pivotal role. An widely used practice to 
extract the temperature of hot nuclear systems is to take 
resort to double-isotope ratio thermometry as suggested 
by Albergo et al [l6j . In an ideal scenario, the primary 
fragments produced in the freeze-out volume are assumed 
to be in their ground states. Particle and 7-decay cor- 
rections to the excited primary fragments have also been 
built in [ItJ. Generally, the feeding effect of secondary 
decay has been accounted through a correction factor 
(15 . HI, [Hi H3] on the measured multiplicities. The tem- 
perature is seen to increase by ~ 10-20 % from the ideal 
situation. 

All these analyzes have been done with the assumption 
that the fragment species produced are noninteracting 
within the freeze-out volume. Strong interaction correc- 
tions, appropriately taken up in the S-matrix approach 
pl| to the grand partition function of the dilute nuclear 
system, where in addition to all the stable mass particles, 
the two-body scattering channels between them can be 
included systematically are seen to modify the fragment 
multiplicities [22|, . The extracted temperature as ob- 



tained in the previous analyses without strong interaction 
corrections may then differ from the real temperature at 
which the fragments were produced. In a schematic cal- 
culation [2~4| in S-matrix approach in dilute infinite nu- 
clear matter, neglecting secondary decay, we found that 
the scattering effects on the extracted temperature and 
volumes are not negligible. In the present communica- 
tion, these ideas are incorporated to provide a realistic 
framework to analyze the data in an experimental multi- 
fragmcntation set up to extract the temperature and vol- 
ume of a finite disassembling nucleus at freeze-out with 
explicit inclusion of 7 and particle decay as well as effects 
from scattering between different fragment species. 

The details of the S-matrix approach, as applied to 
nuclear systems, are given in Refs. [2J, [24|, [25|. For 
completeness, a few relevant equations are presented here 
highlighting the approximations. The grand partition 
function Z of a system in thermodynamic equilibrium 
can be written as a sum of three terms [241 



In Z = In Z gr + In Z° ex + In Z s 



(1) 



The first and second terms correspond to the contribu- 
tions from the ground states and particle-stable excited 
states of all the produced fragment species behaving like 
an ideal quantum gas. The last term sums up the contri- 
butions from the scattering states, expressible in terms 
of the S'-matrix elements. Formal expressions for these 
three terms are spelt out in Ref. [23j . 

The scattering channels, for convenience, can be sepa- 
rated into two parts, one containing only light particles 
and the other the heavy ones, i.e., 



lnZ sc = lnZ' + lnZ* 



(2) 
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The scattering of the heavy ones is dominated by a mul- 
titude of resonances near the threshold; the S'-matrix el- 
ements can then be approximated by resonances, which 



2 



like the excited states are again treated as ideal gas terms 
[26| . These are the particle unstable states. Structurally, 
In Z[} c being then similar to In Z® x , In Z® x and In Z^ c are 
combined together to give lnZ ex (= \nZ® x + lnZ^ c ), 
which contains contributions from particle-unstable ex- 
cited states besides the particle-stable ones. 

In \nZ l sc , only the elastic scattering channels for the 
pairs NN,Nt,N 3 He,Na (N refers to the nucleon) and 
aa have been included. These calculations involve virial 
coefficients [HI, [27j that are functions of only experimen- 
tal entities, namely, phase shifts and binding energies. 
Once the partition function is obtained, total fragment 
multiplicities Yj for the i-th fragment species with N% 
neutrons and Zi protons can be evaluated as 



Yi = Ci 



(3) 



Here Q[= Cz»,JV»] is the effective fugacity defined as 



. B(A il Zi) is the binding en- 



ergy of the fragment and HZi,Nt is its chemical potential, 
which from chemical equilibrium is HZi,Ni = Nifi n +Zifj> P , 
/in and Hp being the neutron and proton chemical poten- 
tials obtained from the conservation of the total neutron 
and proton numbers of the system, (3 is the inverse tem- 
perature. 

For relatively low density and not too low temperature, 
assuming that the quantal distribution can be replaced 
by a classical one, expressions for the primary fragment 
multiplicities of the ith species can be derived as, 

3/2 

y. - V A i Jn n N t +^Zi+B(A t ,Zi)}/T 

1 ~ V A3 



£j=eo 



Y % 



(4) 



In Eq. (4), V is the volume of the system, A = y/2ir/mT 
(we use natural units h — c — 1) is the nucleon ther- 
mal wavelength and g l and g* are the degeneracies of 
the ground and excited states. The sum over the ex- 
cited states includes both 7 and particle-decay (reso- 
nance) channels. In different variants of the models of 
nuclear statistical equilibrium, only the first term (that 
also implicitly contains scattering corrections from res- 
onances in heavy fragments) on the right hand side of 
Eq. (4) has been used to obtain the nuclear thermody- 
namic obscrvablcs. The last term Y s z c is the contribution 
to the fragment yield from scattering, it is nonzero only 
for the fragments in the light species set. Expressions for 
the multiplicity yields Y ™, Y f r etc. (collectively written 
as yj c ) are given in Ref. [22|]. From now on, corrections 
obtained with the use of Y l sc in the extraction of nuclear 
parameters would be called scattering corrections. 

The multiplicities of the primary excited fragments as 
obtained in Eq. (4) undergo changes because of subse- 
quent particle emission. The secondary yield can be writ- 
ten in terms of the variables V, /i„, (J-p and T at freeze-out 



as follows. For light fragments (Ai < 4, Zi < 2), 

3/2 

YJAi, Zi) = vtfJk-^^+z^+B^.z^m + 
A d 



.4 



3/2 

3 p [(N j », n .+Z j », p +B(A j ,Z j ))/T\ 

A 3 



x J* {A j , Zj , T)x? (A, ,Z J ,T)} + Y: c . (5) 

The light fragments are assumed to be produced only in 
their ground states, their multiplicities being given by 
the sum of the first and the last terms in Eq. (5). Their 
population is further fed from decay of heavier species 
given by the second term. The sum j runs over all species 
with Aj > 4 and Zj > 2 having particle-unstable excited 
states and the sum kj runs over all the particle-decaying 

k ■ 

states of the j-th species. The quantity corresponds 
to the branching ratio of the fc^-th state for emitting the 
i-th species; it is calculated using the Weisskopf-Ewing 
model [28| . The quantity Up 3 (Aj , Zj , T) is the internal 
partition function for the particle-unstable states 



■(Aj,Zj,T) = g kj e- e *i /T . 



(6) 



For heavy particles (^4 > 4, Z > 2), the observed yield 



is 



A3/2 



Y(A,Z) = V— e (N^+Zp, p +B(A,Z))/T 

x{g (A)+^(A,Z, T) 



A 



+ a i \3/2 (nifj, n +ZifJ l . D +B(A+ai,Z+z i )-B(A,Z))/T 



(A + a u Z + z h T)x\ 3 [A + a,, Z + z h T)}. (7) 

In Eq. (7), w 7 (= ^2 k gke~ ek ^ T ) is the partition function 
for 7-decaying states, the sum i runs over the emitted 
ejectiles for which we take only n,p, d, i, 3 He and a, Oj, Zi 
being their mass and charge. Kolomiets et al. [17| ar- 
rived also at expressions of the type given in Eqs. (5) 
and (7), the important difference being the absence of 
the scattering correction and consideration of only the 
dominant decay mode. They further considered only nu- 
cleon and a-decay channels. Moreover, the feeding to the 
light fragment yield was neglected. Given a set of exper- 
imental yields for four fragments, their single ratios are 
constructed using Eqs. (5) and (7) resulting in a system 
of three independent equations. The equations are solved 
iteratively in Newton-Raphson method yielding values of 
\x n , Hp and T. The volume can then be determined know- 
ing the yield of a fragment. 

To explore the effect of scattering on the extracted val- 
ues of T and V of a hot fragmenting system, we take 
resort to a numerical experiment. The primary frag- 
ment yields are calculated with given freeze-out tem- 
perature Tf z and volume Vf z in the S'-matrix approach 
as elucidated. The secondary yields are then calcu- 
lated using Eqs. (5) and (7). These are taken as ob- 
served numerical data. In Eq. (7), the first term in the 
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FIG. 1: In panels (a) and (b), the extracted temperature T apP 
for the fragmenting system 124 Sn as a function of the freeze- 
out temperature Tf z shown for two different thermometers 
under different approximations using single ratios (SR). In 
panels (c) and (d) the same are shown using least-squares 
method (LS) for two sets of isotopes as mentioned in the text. 
The dashed-dot, dotted, dashed and full lines correspond to 
Taib, T y , T 1+p and T 1+p+ac = T/~, respectively. 



braces corresponds to the ideal Albergo condition yield- 
ing T a ii,,V a ib, addition of the second term gives the 7- 
decay corrected values T y ,Vy, further addition of the 
last term gives 7 + p (particle)-decay corrected values 
T 1+ p, Vy+p. Only if heavy fragments are taken for mul- 
tiplicity ratios, T 1+p = Tf z and Vy+ P — Vf z . For light 
fragments, to arrive at the actual values of Tf z and Vf z , 
one has to further consider scattering corrections as given 
by the last term in Eq. (5). 

The calculations are done with 124 Sn as a represen- 
tative system. For the fragment species, all the stable 
nuclei upto A=124 and Z=50 as well as their binding 
energies are taken from (29|. All discrete levels upto an 
excitation energy of 20 MeV with lifetimes > 200 fm/c 
as well as their decay modes for 5 < A < 16 have been 
taken [3(| HH into consideration. For still heavier nuclei, 
the sum over excited states in Eqs. (5) and (7) is re- 
placed by an integral convoluted with the single-particle 
level density tu(A, E) [23l. |32|. The integration limits are 
taken between 2 MeV (approximated for the first ex- 
cited state) and 8 MeV (the last particle-stable state) 
for the 7-decaying levels; the resonance limit is taken 
as 20 MeV. The calculations done at different tempera- 
tures in a freeze-out volume 4Vb (Vo is the normal vol- 
ume of 124 Sn) are presented. We have chosen two sets of 
four fragments, namely, 3 He, 4 He, 6 Li, 7 Li and 3 He, 4 He, 
10 Be, 11 Be which we refer to as He-Li and He-Be ther- 
mometers. The extracted apparent temperatures T app 
are found to be quite sensitive to the different approxima- 
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FIG. 2: The same as in Fig. 1 for the extracted volume V app 
in units of freeze-out volume Vf z - 



tions as displayed in the left panels of Fig. 1 . Except for 
Tan,, the other temperatures are not very sensitive to the 
choice of thermometer. Successive improvement of ap- 
proximations is seen to bring the apparent temperature 
closer to the real one. With inclusion of effects due to 
(7 + p)-decay and scattering, the apparent temperature 
T-y+p+sc when calculated yields the actual temperature 
Tfz- The effect of scattering is seen to be substantial. 

The volume V app (measured in units of Vf z ) extracted 
in different approximations is displayed in the left panels 
of Fig. 2 as a function of Tf z for the above mentioned 
thermometers. Scattering has a comparatively more sig- 
nificant role here than that observed in the determination 
of temperature. Its inclusion collapses the apparent vol- 
umes Vapp /Vf z to unity. 

The method so discussed suffers from two limitations. 
For many thermometers, there ma y n ot be convergence 
for the solution as noted earlier [lj]. We also found 
that there may be multiple solutions. We have presented 
those solutions that are robust in the sense that taking a 
considerable range of initial guess values in the iterative 
method, same solutions are obtained. To overcome these 
limitations, we propose that the least-squares fit to the 
secondary multiplicities may be more fruitful in extract- 
ing the temperature and volume. Given experimental 
yields for a chosen number of fragment species n; so , the 
least-squares fit to 



i=l 



Y i {T,V l x n ,n p )] 2 = x 2 



(8) 



has been performed. The quantities are the exper- 
imental multiplicities which arc functions of the thermo- 
dynamic variables at freeze-out and Yi(T, V, fi n , /z p ) are 
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the yields calculated from Eqs. (5) and (7) with various 
approximations as explained earlier. In our calculations, 
Y? xp are taken from our numerical experiment. The ex- 
tracted temperatures T app in the least-squares method for 
the system 124 Sn at different given Tf z and at a freeze- 
out volume Vf z = 4Vo under different approximations are 
displayed in the right panels of Fig. 1. The calculations 
have been performed using a set of light isotopes with 
niso =6 (n,p, d,t, 3 He and 4 He). The calculations are 
repeated with a broader set (rii SO =13) that includes, be- 
sides the light set also the nuclei 6 Li, 7 Li, 9 Be, 10 B, 12 C, 
14 N and 16 O. The 7-decay corrected temperature T 7 is 
found to be insensitive to the choice of fragment set and 
underestimates Tf z (= T 7+p+sc ) considerably. Inclusion 
of particle-decay narrows the gap from Tf z significantly, 
particularly for the broader set of fragment species. Right 
panels of Fig. 2 display the extracted volume V app as a 
function of the freeze-out temperature. The 7-decay cor- 
rected volume Vy overestimates Vf z significantly. Inclu- 
sion of particle-decay brings it closer to Vf z , particularly 
for the larger set. The uncertainty in the (7+p)-corrected 
value for the volume, with n.; s . =13, is seen to be at most 
25% and that for temperature, it is at most 5%. Inclu- 
sion of heavier species in the fitting procedure masks the 
scattering effects. The calculations have been repeated 
for Vf z =SVo and 8Vb; the conclusions do not change for 
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